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ABSTRACT 

We study the spatial distribution of satellite galaxies around isolated primaries using 
the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. 
We select isolated primaries from the spectroscopic sample and search for potential 
satellites in the much deeper photometric sample. For specific luminosity primaries 
we obtain robust statistical results by stacking as many as 50, 000 galaxy systems. 
We find no evidence for any anistropy in the satellite galaxy distribution relative 
to the major axes of the primaries. We derive accurate projected number density 
profiles of satellites down to 4 magnitudes fainter than their primaries. We find the 
normalized satellite profiles generally have a universal form and can be well fitted by 
projected NFW profiles. The NEW concentration parameter increases with decreasing 
satellite luminosity while being independent of the luminosity of the primary except 
for very bright primaries. The profiles of the faintest satellites show deviations from the 
NFW form with an excess at small galactocentric projected distances. In addition, we 
quantify how the radial distribution of satellites depends on the colour of the satellites 
and on the colour and concentration of their primaries. 

Key words: Galaxies: dwarf, Galaxies: structure, Galaxies: Local Group, Galaxies: 
fundamental parameters 



1 INTRODUCTION 

In the ACDM model, smaller structures falling into larger 
haloes can survive there as substructures and host observed 
satellite galaxies. Therefore, the distribution of satellites 
around primaries holds important information about galaxy 
formation, the population of substructures and even the 
nature of dark matter. In the past decade or so, fainter 
satellites around the Milky Way (MW) and M31 have 
been discovered in the Sloan Die 
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The large body of work on satellite galaxies reflects the 
fact that they are not only a critical small scale test of the 
ACDM model, but also a probe of the nature of dark matter; 
yet the satellite data with which all theories and models are 
compared with are merely those of two primaries, the MW 
and M31. Although the satellite populations of MW and 
M31 are known better than other satellite systems, there 
is no guarantee that they are typical. Clearly, robust and 
reliable tests of cosmological and galaxy formation models 
require comparison with a statistically representative sample 
of galaxies and their satellites. 

Early studies were limited by the relatively small 
satellite samples av ailable at the time Hlo lmbcrg 196^; 
iLorrimer et al.Hl99l ; [Zaritskv et al.lll993l . [l997bl '). With the 
advent of large galaxy redshift s urveys such as the 2dF 
Galaxy Red shift Survey (2dF GRS: IColless et"al]|200lD and 
the SDSS l|York et all |2000| ). it is now possible to con- 
struct external galaxy samples spanning a much larger vol- 
ume. Studies with signiflcantly improved statistics have 
been carried ou t usin g these new surveys. For example, 
ISales fc LambasI (|2004l ) studied the spatial distribution of 
satellites ar o und p rimaries using the 2dFGRS. More recently 
lYang et al. I (|2006l ) studied how spectroscopically identified 
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satellite galaxies were distributed in SDSS groups relative 
to the orientation of the central galaxy. However, due to the 
flux limit of redshift surveys, analyzing the satellite systems 
of external isolated galaxies is still challenging because typ- 
ically only one or two satellites are detected per primary 
galaxy. In addition, the real space position of satellites with 
respect to their primaries is uncertain due to redshift space 
distortions and projection effects. To circu mvent the afore- 
mentioned problems, we ||Guo et aLlboill . hereafter Paper 
I) have developed a method of stacking the primaries and 
their satellites in order to obtain a fair and complete sample 
that can yield statistically robust results for certain classes 
of primary galaxies. This method has been successfully ap- 
plied to the estimation of the satellite luminosity functions 
of isolated primary galaxies in the SDSS. 

In this work, for the same primary and satellite sam- 
ples we explored in Paper I, we are now interested in the 
average spatial profile of the distribution of these satellites 
around their primaries. These density profiles are an im- 
portant tracer of the distribution of substructures in the 
primary halo and can provide us with useful information to 
test current models of the formation and evolution of dark 
matter haloes. 

In the cold dark matter (CDM) cosmological model, 
the density profiles of dark matter haloes follo w a uni- 
versal form ijNavarro. Frenk. fc Wh"it3l 19961 . 1 19971 . hereafter 
NFW profiles) with an inner cusp, p(r) oc r~^, and an 
outer slope of p(r) oc . The transition scale, rg, is nor- 
mally specified through the concentration, c = r2oo/''s, 
where r2oo is defined as the radius enclosing a mean interior 
density 200 times the critical density. Besides the overall 
mass profile, it is remarkable that the spatial distribution 
of dark matter substructures, which could host satellites 
galaxies, also follows this universal form independent of the 
mass of the substructures (iDiemand. Moore, k, Stadelll2004l : 
ISpringel et ai]|2008l : iLudlow et al.ll2009l '). However, the num- 
ber of observed satellite galaxies around the MW and M31 
is much sma ller than th e numbe r of substruct ures predicted 
by ACDM (iKlvpin et al. 1999: ' Moore et akl figgg). giving 
rise to the so-called "missing satellites problem". Statisti- 
cally robust number density profiles of observed satellites 
will certainly help us understand how satellite galaxies pop- 
ulate the substructures. In addition, a reliable density profile 
is required to extrapolate the incomplete observational data 
of satellites around the MW and M31 to compare with mod- 



els (e.g. iKoposov et"al]|2008l : [Tbllerud et al.ll2008l ). 

The recognition of the importance of the spatial profiles 
of systems of satellite galaxies has resulted in many studies. 
Early work with samples from a limited volume have mainly 
focused on fitting the slope of the densi ty profile of satel- 
lite galaxies a round isolated primaries ("Lake fc Tremaing 
Il98d : IVader fc Sa ndagc 1991; Lorrimcr ct al. 1994). With 
large galaxy redshift surveys, the dependence of the pro- 
files on the colour and morphology of prima r ies has be- 
gun t o be explored (e.g . ISales fc Lambajr2005l : IChen et "al] 
l2006l 'l. iKlvpin fc Pradal l|2009l l studied the projected num- 
ber density profiles and velocity dispersion aro und isolated 
red pr imaries using the SDSS redshift sample. iMore et al.l 
l|2009h used an iterative method, tested on mock galaxy 
catalogues, to find satellite systems around central galaxies 
with a range of luminosities in the SDSS. The distribution 
of velocities of the satellites was used to infer mass-to-light 



ratios as a function of central galaxy luminosity. Closely re- 
lated to this are studies of t he radia l distributions of satellite 
galaxies in clusters, gro ups |Li et al. 20 07: Wang ct al. 20li]) 
and on smaller scales jWatson et al.l I2OI2I '). Further work 
has focused on the distri bution of satellites aro und inter- 
mediate re dshift galaxies (iNierenberg et al.ll201ll ). elliptical 
primaries jSmith. Martinez, fc Graham 20041') and isolate d 
galaxies in the SDSS ( Lares. Lambas. fc Dominguej[201ll ). 
Besides the stud ies that sta ti stical ly estimate mean number 
density profiles, iKim et al] (|201ll ) have directly measured 
the number density profile of the nearby field galaxy M106. 

In this paper, we are not only interested in the projected 
number density profile of satellites around isolated primary 
galaxies binned by luminosity, colour and morphology, but 
also on the dependence of the profiles on the properties of the 
satellites themselves. To this end, we select our primary sam- 
ples from the SDSS DR8 spectroscopic sample (~ 660 000 
galaxies) and satellite candidates from the photometric sam- 
ples (~ 96 000 000 galaxies) with the same criteria as in 
Paper I to build significantly large samples. We restrict the 
photometric sample to galaxies brighter than rrir = 20.5 as 
in Paper I (See Paper I, Section 4) to ensure completeness. 
Based on these large samples, we explore the dependence 
of the density profiles on the properties of primaries and 
satellites. 

The remainder of this paper is organised as follows. In 
Section 2, we briefiy describe the selection of primary galax- 
ies and their satellites; in Section 3, we develop a method 
for estimating the projected satellite number density pro- 
file; in Section 4, we present our estimate of the projected 
satellite number density for different primary samples. We 
conclude, in Section 5, with a summary and discussion of 
our results. Throughout the paper we assume a fiducial 
ACDM cosmological model with Qm = 0.3, Qa = 0.7 and 
Ho = 70km s~^ Mpc-\ 



2 SAMPLE AND METHOD 
2.1 Data and sample selection 

We have built two separate catalogues similar to those in 
Paper I. The smaller one is of galaxies from the main SDSS 
spectroscopic catalogue from which we select our primary 
galaxies (hereafter the spectroscopic catalogue). The larger 
one is of galaxies with photometric redshifts and magnitudes 
from which we select the neighbouring galaxies (hereafter 
the photometric catalogue). The spectroscopic catalogue is 
constructed from the SDSS DR8 spectroscopic subsample 
(north galactic cap) including all objects with high qual- 
ity redshifts (zconf > 0.7 and specClass = 2) and a Pet- 
rosian magnitude r ^ 17.77. The photometric galaxy cata- 
logue is from the SDSS DR8 photometric subsample (north 
galactic cap) and includes only objects that have photomet- 
ric redshifts, none of the flags BRIGHT, SATURATED, or 
SATUR_CENTEE0 set and model magnitudes r sC 22.0. We 



^ When applying our isolation criteria to reject primaries with 
bright neighbours we use a source catalogue that also includes 
objects for which SATURATED and/or SATUR.CENTER flags 
are set. These objects are mainly stars and we prefer to reject 
systems contaminated by bright stars as the presence of such 



© 0000 RAS, MNRAS 000, 000-000 



Satellite Number Density Profiles 3 




Figure 1. Schematic diagram of the sample selection procedure. 
For each acceptable primary, there should be no neighbouring, 
sufficiently bright, galaxies within a cylinder of radius -Riso, cen- 
tred on the primary, and nearby in redshift space. As defined in 
Fig. 1 in Paper I, nearby means either \zc — 2:" '^'^''1 < Azg, where 
Azb is the maximum allowed spectroscopic redshift difference be- 
tween a primary (c) and another galaxy (s), or l^c — Zp'^'^'^l < 
apCTp, where tTp is the measurement error of the photometric red- 
shift and Op is the tolerance of the error, for galaxies that have no 
spectroscopic redshift. Satellites will lie nearby in redshift space 
and within the cylinder of radius -Rinnor (red), whereas the lo- 
cal background to be subtracted is determined from the volume 
between this inner cylinder and the outer one with radius -Router ■ 

select only objects with corresponding entries in the SDSS 
database PhotoZ table, which naturally selects galaxies and 
excludes stars. As galaxies with r ^ 17.77 are included in 
both SDSS catalogues, a small fraction of the photomet- 
ric catalogue galaxies also have spectroscopic redshifts. We 
use de-reddened model ugriz magnitude s and fc-correct all 
galaxi es to z = Q with the IDL code of iBlanton fc Roweij 
|2007l ). We estimate l/-band magnitudes from q and r-band 
magn itudes assuming V = gi — 0.55((7 — r) — 0.03 (|Smith et al.l 
I2OO2I ) and all our sample selection and magnitude cuts are 
performed using this l/-band magnitude. 

Our sample of isolated primary galaxies is chosen us- 
ing the same criteria as in Paper I, illustrated in Fig. [1] 
We select primary galaxy candidates of absolute magnitude, 
Mp, in the range Mc - AMbin < Mp < Mc -I- AMbin. 
We then filter these primary candidates, using a series of 
criteria summarised in figure 1 of Paper I, to guarantee 
that a) there are no luminous neighbouring galaxies pro- 
jected within 27?inner of the primary, unless these luminous 
neighbours are sufficiently separated in redshift from the 
primary and appear here due to a chance projection; b) 
the satellite search areas (projected distance i?inner from 
the primary) around each primary do not overlap with 
each other. Further details of the generation of the two 
samples are may be found in Paper I. The values of the 
selection parameters, {AA/bin, AA/faint, Azs, Qp, jti^'™} = 

unmasked stars could effect the efficiency with which background 
galaxies are detected. 



{0.5,0.5,0.002,2.5,20.5}, are the same as the default val- 
ues in Paper I. Here AAfbin is the half-width of the primary 
magnitude bin, AAffaint is the magnitude difference be- 
tween the primaries and satellites used to isolate primaries, 
Azs and Qp are the parameters used to exclude galaxies 
that are at a significantly different spectroscopic redshift 
and photometric redshift respectively. The meaning of these 
parameters is explained in Fig. [1] (see section 2 of Paper I 
for more details). One small change relative to Paper 1 is 
that -Rinnor is choscu to increase with increasing primary lu- 
minosity, in order to ensure that no satellites are missed for 
the most luminous primaries. 

2.2 Estimating the projected satellite number 
density profile 

Once our primary galaxies are defined, their potential satel- 
lites are found from the photometric galaxy catalogue as 
depicted in Fig. [2] The method we use is similar to that in 
Paper I, except that for each primary, the number of galax- 
ies is counted and binned by their projected radial distances 
from the primary, -/V(r|""), as well as by their luminosity. 
That is, for the ith primary galaxy, the number of inner 
galaxies in each annulus, A''j'""°'^(r|""), is found by counting 
all the neighbouring galaxies within the annulus of radius 
r™'' that satisfy the following conditions: at least AMfaint 
fainter than the primary; if they have a spectroscopic red- 
shift, Zs, then they should satisfy |2c — ^s] < Azs; or if they 
only have a photometric redshift z^, then they should satisfy 
\zc — Zp\ < QpCTp, where cr* is the error in the photometric 
redshift as defined in Paper 1. The number of outer galax- 
ies, Af""*"^"^, is determined by applying the same conditions 
to galaxies in the outer area between 7?inner and -Router. As- 
suming, for now, that few genuine satellites will be projected 
beyond -Rinnor we can estimate the surface density of genuine 
satellites in each annulus as 

jyinner/ ann\ ivmuter 
j-jsat/ ann-, _ i \' j J _ i Q^ 
^ ^ J ' ^ann pouter ' V ' 

where j4f"" and ^4°"''^'' are the areas of the inner annulus and 
outer region respectively. If necessary, we take account of the 
sky coverage of SDSS DR8 by reducing the areas and 
p outer ^^^g amounts defined by the DR7 mask described 
in iNorberg et all (|201ll ). 

Because of the apparent magnitude limit of SDSS, 
which we take to be m'™ — 20.5, we count the faintest 
satellites only at low redshift and are progressively limited 
to more and more luminous satellites with increasing red- 
shift. To account for this and construct an unbiased estimate 
of the projected satellite number density profile over all pri- 
mary galaxies, we accumulate the area contributed by the 
ith primary to the j'th annulus for the detection of satellites 
brighter than A^slt"", 

r A - A//"'"^"" ^ A/7"'™ 

/ ..rtrun-, _ J ^i] JWsat ^ ■'"i /r,\ 

A,,(Ai,,t ^ - I M',ir > Mf^ . ^ > 

Here Aij is the unmasked area of the jth annulus surround- 
ing the ith primary and Af|™ is the absolute magnitude 
that corresponds to the apparent magnitude limit, m'™, 
of the SDSS photometric catalogue at the redshift of the 
primary. A/]™ = m'™ - 51ogio(A^) - K{zi), where 
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Figure 2. Schematic sliowing tiie selection of potential satellite galaxies in annuli of radii r^"" within -Rinncr a-^d a- reference sample 
bounded by the radii -Rinner < r < -Router, used to subtract the residual contaminating background. For both samples we apply the 
stated redshift cuts to reduce background contamination. We also apply the stated absolute magnitude cut to both samples (assuming 
the neighbouring galaxies are at the same redshift as the primary) though this cut is redundant unless Router > SRinnor as otherwise 
the existence of such bright neighbouring galaxies would automatically lead to the exclusion of the primary galaxy. 



and K(zi) are the corresponding luminosity distance and 
^-correction. This contributing area is set to zero if any po- 
tential satellites within the magnitude bin are too faint to 
be included, in which case we exclude this primary and its 
satellites as its contribution to the mean projected satel- 
lite number density profile would be incomplete. We fur- 
ther define NffiM^lt") to be the number of detected po- 
tential satellites brighter than Afga"" in the jth annulus sur- 
rounding the ith primary and A^i^*(Msat") to be the corre- 
sponding number of detected galaxies in the outer annulus, 
J^innor < r < Router, whosc uumasked area is Hence 
we can express the mean surface density of satellite galaxies 
brighter than A/g'^"" in the jth annulus as 



j..at(^arr„^^truu) 



(3) 



In practice, we divide the projected radial distance from 
the primary into 20 bins (j = 1, 2, • • • , 20). Because of a con- 
cern that the SDSS data reduction pipeline may occasionally 
misclassifies fragments of the spiral arms of bright galaxies 
as separate galaxies we exclude individual annuli that are 
within 1.5 times the Petrosian radius, -R90, of the primary 
galaxy. We set our magnitude limit, Mga"", either by abso- 
lute value, such as —20, —19, —18, or by magnitude relative 
to the corresponding primary, M*^t"" = Mp + 1.0, 2.0, 3.0. 
We present results using both thresholds so that we can de- 
termine whether the number density profiles depend on the 
absolute luminosity of satellites or on the relative luminosity 
between satellites and their primaries. 

The process of estimating the projected satellite num- 
ber density profiles is quite similar to estimating the 
satellite luminosity functions in Paper I. We divide our 
primaries into three luminosity bins centred on Mc = 
—21.25, —22.0, —23.0. The choice of the parameters -Rinncr 
and -Router are a balance between making them sufficiently 
large to avoid severely truncating the density profiles and 
making them too large such that our sample shrinks due to 
the selection process excluding overlapping systems. A sen- 



sible choice is to set -Rinner to exceed the anticipated size 
of the satellite system, r20(£|, and -Router to be roughly a 
factor of two larger so as to get a good, but still local, es- 
timate of the background density. Here we have adopted 

the following values of (-Rinner, i?outer), (0.3, 0.6) , (0.4, 0.8) 

and (0.55, 0.9) Mpc for primaries in magnitude bins, Ale = 
—21.25, —22.0 and —23.0, respectively. The values of -Rinner 
have been compared with the mean of the estimated r2oo 
values for each galaxy in the chosen magnitude bin (see Sec- 
tion 3) and are found all to be larger, suggesting that the 
search radii for the different primary magnitude bins are 
sufficiently large to capture all satellites. Additional reas- 
surance is provided by the tests in Appendix [C] which show 
that the profiles are insensitive to changes in the values of 

dinner Or -Router. 



2.3 Exploring the angular distribution of satellites 

The projected radial satellite number density profile, 
E^^'(r), which is the focus of this paper, is the azimuthal 
average of the 2D surface density, 'E^'^^{r,6), where 6 can 
be taken as the position angle between the major axis of 
the primary and the line connecting the primary and satel- 
lite (see Fig [2}. The angular dependence of this distribu- 
tion may also carry information on the formation and evo- 
lution of the satellites around their primaries. For exam- 
ple, if we assume that the satellite galaxies inhabit an un- 
biased set of dark matter subhalos, then we would expect 
satellites to cluster preferentia lly along the major axis of 
the halo (|Libeskind et al.ll2005f ). Moreover, it is known that 
the host halos of satellite systems are accreted from fila- 
ments, which can cause the angular distribution of satellites 
around primaries to be anisotropic (|Hartwickll2000l ). In fact, 
numerous such anisotropies have been observed. For exam- 
ple, the famous "Holmberg Effect" (jHolmberd Il969l ') . sug- 



^ Here r2oo depicts the radius at which the mean interior density 
is 200 times the cosmological critical density. 
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Figure 3. Schematic sliowing the definition of the position angle 
6 which characterises the angular position of the satellite relative 
to the direction of the primary galaxy's major axis. The grey and 
blue ellipses are the primary and satellite galaxy, respectively. 
The angle ipp is the position angle of the primary. 



gests satellites of isolated, large, and inclined spiral galax- 
ies are preferentially located along the minor axe s of thei r 
primarie s, a result s uppor t ed by IZaritsky et al] (Il997bl). 
Howev er 'Yang et al.' (2006') , 'Azzaro et all (|2007^ ■ iBrainerdI 
l|2005l ) and Agustsson & Braincrd ( 201ll) found the opposite 
effect that satellites prefer alignment with the major axis, 
especially for the satellites of red primaries. Since the direct 
observation of satellite systems is not easy, the sample of ex- 
ternal satellite systems is limited in both volume and quality. 
These contradictory results may suggest that the mean am- 
plitude of the anisotropy could be very weak, or the form of 
the anisotropy could be dependent on the selection of the 
primaries or even the satellites themselves. 

It is therefore interesting to attempt to quantify the 
mean anisotropy of our large sample of satellite galaxies. We 
characterise the angular distribution using the position an- 
gle 9 described in Fig. O The assumed elliptical symmetry of 
the primary implies that 9 ranges from 0° to 90° with these 
extremes indicating that the satellites are located along the 
major or minor axis respectively. The anisotropy of the an- 
gular distribution is then quantified by the probability dis- 
tribution of the angle 9. In practice, accurate measurement 
of 9 requires a robust measurement of the position angle (f)p 
defining the orientation of the prim aries. To achieve this, we 
adopt the same selection criteria as lSiverd. Rvden. fc GaudU 
l|2009D . We only select primaries and their satellites that sat- 
isfy the condition, qiso < 0.9 and gmom < 0.9, where giso is 
the isophotal axis ratio defined as giso ~ aiso/6iso and Qmom is 
the adaptive moments a xis ratio, gmo m = ((1 — e)/(l-l-e))^'^^ 
where e = {el + el fl'^ l|Rvdenll200-l FI. We also exclude the 
primaries together with their satellites if there is a discrep- 
ancy of more than 15 degrees between the measured isopho- 



tal and de Vaucouleurs position angles, A6p 



> 15° 



and ex are second-order parameters from SDSS, where r : 



M^3;-|-Myj;,e+ = {M^x- Myy)/T, 

are the second-order adaptive moments. 




Figure 4. The probability distribution, P{6), of the position an- 
gle 6 of satellite galaxies. In the top panel the blue and red lines 
show the distributions, P""'°''((?), of inner galaxies, and P°"*™((?), 
of outer galaxies, respectively. The inferred distribution for true 
satellite galaxies, P'^'^*{9), is shown by the black line. The lower 
panel shows the distributions for red and blue primary subsam- 
ples, as the red and blue solid lines respectively. The error bars 
are calculated from 1000 bootstrap resamplings. The expectation 
for a uniform distribution is shown by the dashed lines in the two 
panels. 



For our sample of selected satellite systems, the number 
of satellites located at angle 9j around the ith primary, can 
be estimated as 



Nr\9j) = Nl'''"''(9,) 



A° 



(ivr^^(e,)) 



(4) 



where (A^°"''^'^(^j)) represents the azimuthal average as we 
assume the background galaxies are, on average, isotrop- 
ically distributed. Note that to avoid biasing the angular 
distribution we exclude systems that are incomplete due to 
the survey mask. We can then define an unbiased estimator 
of the average distribution of the satellite position angle 9 
for all selected primaries as 



N{9,,Ml 



) 



(5) 



■ 2Mxy/T and Mxx here 



The normalized probability distribution of 9 for pri- 
maries in the magnitude bin Mc = —22.0 is shown in 
Fig. [l] The blue, red and black solid lines in the top panel 
are the probability distributions of 9 for inner galaxies, 
P'""'='(6l), outer galaxies, P°"*"(e), and that inferred for 
satellite galaxies, P{9), respectively. In the bottom panel of 
Fig.[4l we show the probability distribution of 9 for two satel- 
lite subsamples split by rest frame colour of their primary 
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Figure 5. The mean projected number density profiles of satellites around primaries for various selections of primaries and satellites. 
The left hand panels (a and b) show the mean number densities and the right hand panels (c and d) show the same profiles but now 
normalized by the total number of satellites within the r2oo radius (see right-hand axis) and with the radius expressed in units of the 
adopted mean r200 radius of the primaries. In each panel the different coloured lines correspond to primaries of differing luminosity as 
indicated in the legend. The profiles in the upper panels (a and c) are for satellite samples brighter than M*^"° = —19, while the lower 
panels (b and d) are for satellites that are less than 2.0 magnitudes fainter than their respective primaries. The error bars show the errors 
on the mean profiles estimated by bootstrap resampling. 



galaxies. These samples are divided according to the well- 
known colour bimodality in the colour- magnitude plane (e.g . 
IStratcva ct al."200l': 'Baldrv et aLll2004l : IZehavi et al]|2005l ). 
Following Zchavi ct al. ( 20ol)7 we use an equivalent colour 
criterion of ° °(q - r)cut = 0.19 - 0.24Mr (not identical to 
IZehavi et al.l as our magnitudes are fc-corrected to z = 0.0 
rather than z = 0.1). The probability distributions of 6 in 
Fig. [4] are all consistent with isotropic distributions. This is 
confirmed by a two-sample Kolmogorov-Smirnov (KS) test 
that compares the distributions of inner and outer galax- 
ies for the whole primary sample and also the subsamples 
of red and blue primaries. The two-sample KS probabilities 
from the whole sample, blue and red primary subsamples 
are 0.64, 0.37, 0.87 respectively, which implies that the pairs 
of distributions have no statistically significant differences. 
The same tests for primaries in other magnitude bins show 
similar results. 



Therefore, with our satellite system sample, we find that 
there is no statistically significant evidence that the distribu- 
tion of satellites around primaries is anisotropic. This could 
signify that the anisotropy of the distribution of satellites 
around isolated primaries is intrinsically insignificant. How- 
ever, one also has to keep in mind that our inner sample 
includes contamination from interlopers, because these are 
only rejected using inaccurate photometric redshifts, and 
this could dilute any intrinsic anisotropy signal. 



3 RESULTS 

We now return to the azimuthally averaged density profiles. 
Density profiles for satellites brighter than Mg^t " around pri- 
maries of magnitude Mc are shown in Fig. [S] for a variety of 
primary magnitude bins and satellite magnitude cuts. Pan- 
els [5^ and [Sja show that the number of satellites increases 
with increasing primary luminosity and extends to larger 
radii. To investigate the variation in profile shape between 
different subsamples of satellites and primaries, it is helpful 
to use scaled variables. To this end, we recast the profiles in 
terms of a; = r jr^oo and divide the number densities by the 
total number of satellites within r2oo- 

The values of r2oo used to scale the radii can be de- 
termined from the stellar masses, themselves inferred from 
the measured galaxy luminos ities and co l ours, and the abun- 
dance matching technique of lGuo et al.l (|2010l ). which gives 




where c = 0.129, Mo = lO" * M0, a = 0.926, P = 0.261 
and 7 = 2.440 are fitted constants. The halo mass can be 
related to a radius through 



Mhalo = -g-200pcritr^oo- (7) 

For the primary galaxies there is a significant uncer- 
tainty in the stellar mass that is inferred from the measured 
luminosities and colours. Thus, rather than using the in- 
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Figure 6. The dependence of the scaled satellite density profiles on satellite luminosity for primaries in the magnitude range —22.5 < 
Mp < —21.5. The different panels show different satellite selections. The upper panels, which are very similar, show profiles for satellites 
brighter than a threshold that is either a fixed value (panel a) or specified as a magnitude difference with respect to the corresponding 
primary (panel c; see legend). The lower two panels show profiles for satellites in bands of magnitude again either specified between fixed 
values (panel b) or between values relative to the corresponding primary (panel d). 



dividual r2oo values to normalise satellite number density 
profiles for each primary before stacking them together, the 
mean r2oo for all primaries in the luminosity bin of interest 
is determined and a single rescaling is performed on the un- 
sealed, stacked profile. Using r2oo values determined in this 
way, the = —21.25 and —22.0 samples line up very well, 
as shown in panels [5j; and[5jl. However, the Mc — —23.0 re- 
sults, not shown in Fig. O are slightly offset. For these, the 
most luminous primaries, the relation between stellar mass 
and halo mass becomes very fiat and there is a large spread 
in the halo mass corresponding to a given stellar mass. This 
makes the assignment of a value of r2oo to these primaries 
extremely uncertain. The directly inferred virial radius is 
0.73 Mpc but we find that a smaller value of 0.52 Mpc results 
in better scaUngs. Given the large uncertainty in this assign- 
ment, it is not unreasonable to adopt this smaller value. We 
shall do this in what follows but this uncertainty must be 
borne in mind when interpreting the results for the bright- 
est primary bin. The final values of r2oo for the three pri- 
mary magnitude bins are 0.24, 0.37, 0.52 Mpc, the first two 
of which come directly from equation ((7|. By adopting these 
mean r2oo values, we have a mass-to-light ratio increasing 
with luminosity as M/Ly oc Ly^^ . This is similar to the 
relation found bv lPrada et al.l (|2003l '). albeit for the B-band, 
from a set of spectroscopically selected satellites from SDSS. 
For Mv =-22, we find that M/Lv' = 140. 

First, we explore the dependence of the normalized 
profiles on the luminosity of the satellites. In Fig. [6l we 



show normalized profiles for primaries of fixed luminosity 
(—22.5 < A/p < —21.5) with a variety of different satel- 
lite selections. In all cases we find the outer shapes of the 
density profiles to be very similar. The only variation is on 
small scales (roughly r/r2oo < 0.1) where the density profile 
is steeper and higher for the faintest satellites. 

We examine the shape of these density profiles 
more systematically in Fig. [7] where we fit the un- 
normalized density profiles using an analytic model. 
We have chosen to fit our satellite profiles using 
NFW profiles as they a re known to be good fits to 
both dark matter haloes (jNavarro. Frenk. fc Whitd Il996l . 
[l993) and to the distribution of substructures within 



them ( Diem and. M oore, fc StadeJ2004l : ISpringel et al ]|2008l : 
lLudlow"et al. 2009). To perform these fits we first project 
the NFW profile and subtract from it the mean density in 
an outer "background" annulus as described in Appendix 
[Al so as to match the background subtraction process that 
we have applied to our observational data. We then per- 
form a fit using E'"'*(r) = At,{r,c,r2oo) / M200, where A is 
a scale factor, c is the concentration, and E(r, c, r2oo) is the 
projected NFW profile with background subtracted as given 
in Eqn. IA5I The r2oo radius is fixed to the respective val- 
ues that we have adopted for each of our bins of primary 
magnitude. 

Fig.[7]shows the resulting fits. The bright satellites, such 
as satellites with magnitudes —21.0 < Ms < —20.0 or Mp -|- 
1.0 < Ms < Mp + 2.0 are very well fitted by the projected 
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Figure 7. Fits to the satellite density profiles of primaries of magnitude —22.5 < Mp < —21.5 using projected, background-subtracted 
NFW profiles. The various panels show different selections of satellites as in Fig. [G] except that here we have not normalized the profiles, 
but instead allowed the amplitude, A, of the fitted profiles to float. The measured profiles are shown by the data points and the best 
fitting NFW profiles arc plotted as solid lines. The best fitting amplitudes and concentrations are listed in the legends in each panel. 



and background subtracted NFW profiles. The NFW fits 
also remain good descriptions of the data for the cumulative 
samples of satellites defined by a faint magnitude threshold. 
For these samples, shown in the upper panels of Fig. [T] the 
concentration increases steadily with decreasing luminosity. 
In the lower panels of Fig. [71 which show density profiles 
for satellites in differential bands of luminosity, we see both 
a stronger dependence of concentration on luminosity and 
small deviations from the NFW form for the faintest satellite 
samples. 

We now turn to the dependence of the satellite profiles 
on the magnitude of the primaries. In Fig. [5] we show fits of 
the projected and background-subtracted NFW profiles to 
satellite profiles of primaries in each of our three magnitude 
bins. Each of the panels corresponds to a different selection 
of satellites. We see that NFW fits are good descriptions of 
the satellite distribution regardless of the luminosity of the 
primary. The right hand panels of Fig. [5] show the density 
profiles and fits for sets of satellites defined by fixed offsets in 
magnitude from the magnitude of their respective primary. 
If the combined primary and satellite systems scaled in a 
self-similar way we would expect the three density profiles 
in each of these panels to lie on top of each other. In contrast, 
in each panel, we see systematic variations in the shape and 
amplitude of the profiles with the primary luminosity. If in- 
stead we look at the left hand panels, which show satellites 
selected in different fixed magnitude bands, then we see that 
the concentration decreases steeply with increasing satellite 
luminosity, but is less dependent on the luminosity of the 



primary. With the exception of the brightest primary mag- 
nitude bin (Mc = —23.0), the satellites of a given luminosity 
are more or less distributed in the same way about primaries 
of different luminosity. Only the normalization of the profile, 
as parameterized by A, increases with increasing primary lu- 
minosity. This depends quite strongly on luminosity going 
roughly as the luminosity to the power of 2.5. If we normal- 
ized each of these satellite profiles as we did in Fig. [S] then 
their shapes would show very little variation with primary 
luminosity. 

For the case of primaries in the Mc = —23.0 bin, the 
measured concentrations for the satellite distributions are 
systematically lower than those of similar luminosity satel- 
lites around less luminous primaries. The reason for this is 
not clear. It may be that this result reflects the actual satel- 
lite distribution around bright primaries. However, it could 
be an artefact of the estimation procedure. One possibility is 
that the non-line arity of the halo mass-stellar mass relation 
dGuo et al.ll2010l ) means that the range of actual halo mass 
increases in the brightest primary bin. Thus, stacking all pri- 
maries using a single r2oo value may be introducing errors 
that would smear out the resulting profile. Another potential 
source of systematic error comes from the tendency for faint 
satellites to be missed around bri ght primaries because of in - 
accuracies in the sky-subtraction l|Mandelbaum et al.ll2005l ') . 
Even with the updated sky-subtraction algorithm employed 
for DR8 there may still be some residual loss of faint satel- 
lites around the brightest primary galaxies l|Aihara et al.l 
[20i3). 
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Figure 8. NFW fits to the density profiles of satellites around primaries of different luminosity. Each panel shows fits for the three 
different bins of primary luminosity indicated in the legends. The different panels correspond to different selections of satellites. The 
panels on the left, a), b) and c), are for successively fainter bands of satellite luminosity as indicated on the legend, with the range of 
satellite magnitudes being the same for each primary. In contrast the panels on the right, d), e) and f), are for bins of satellite magnitude 
that are specified as an offset relative to the magnitude of their respective primary. 



3.1 Colour and type dependence 

Our large sample of satellite systems enables us to divide our 
samples by the colour or the type of the primaries. Fig. [9] 
shows the resulting profiles when primaries of l/-band mag- 
nitude —22.0 ±0.5 are split by colour and by concentration. 
Panels a) and c) of Fig. |9]show that the normalized profile of 
satellites around blue primaries is more concentrated than 
that around red primaries. Panels b) and d) split the sam- 
ple into early and late types, where early type is defined as 
having a concentration index C = petroRgo / petroRz,o ^ 2.6, 
with petroRgo and petroR^o being the SDSS Petrosian 90% 
and 50% hght radii respectively. This division roughly sepa- 
rates early type (E/SO) from late-type (Sa/b/c, Irr) galax;ies 
jShimasaku et aLlbOOlh . We also see the amplitude of the 
profiles of late types is suppressed with respect to that of 
the early types. However, the concentration indices, c, from 



the fits similarly show that the concentration of satellites 
around late types is higher than that of early types. 

We can also use the colour information available in 
SDSS to probe the properties of the satellites. Firstly, for 
the bin of primary magnitude, A/c = —22.0, we divide the 
satellites into two luminosity bins, —21.0 < Ms < —20.0 and 
—20.0 < Ms < — 19 and into red and blue subsamples using 
the same cut as before. Fig. llOb .b.d and e show the measured 
profiles of these blue and red satellites and NFW fits. We 
first note from Fig. IIOH and e that for these relatively bright 
satellite samples the abundance of blue satellites is greater 
than that of red satellites at all radii, with this difference in- 
creasing for the fainter sample. The profiles of the brighter 
satellites have a similar shape for red and blue satellites, 
while the fainter red satellites have an excess at ~ 0.1r2oo 
relative to the fainter blue satellites. To investigate whether 
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Figure 9. The satellite profiles for primary galaxies of magnitude —22.5 < A/p < —21.5 split by the type (concentration) and colour of 
the primary. Panels a) and b) show the normalized profiles while panels c) and d) show NFW fits to the unnormalized profiles (see right 
hand axis). In the upper panels the blue lines refer to blue primaries and red lines to red primaries, while in the lower panels blue refers 
to late-type primaries and red to early types. 



these differences are driven by the colours of the associated 
primary galaxies we further split the satellites brighter than 
j^^trun ^ _2Q ^jjg colour of their primary. The results 
are shown in Fig. llOfc and f. Both red and blue primaries 
have more blue than red satellites. The concentrations of 
red and blue satellites around blue primaries are similar. 
Red primaries have lower concentrations for both their red 
and blue satellites, with the blue satellites having a particu- 
larly low concentration. The colour of the primary appears 
to be more important than that of the satellite in determin- 
ing the concentration of the satellites. As shown in Fig. |8l 
the satellite luminosity also has a strong effect. 



4 DISCUSSION 

Using a stacking analysis we have estimated the mean pro- 
jected density profiles of satellite galaxies around a large 
sample of isolated primary galaxies selected from the SDSS 
DR8 spectroscopic galaxy catalogue and we have quantified 
how they depend on the properties of the satellites and pri- 
maries. The selection of primaries and the local background 
subtraction technique, which ma kes use of photom etric red- 
shifts, is the same as in Paper I (|Guo et al.ll201ll ) where we 
estimated the mean satellite luminosity functions of these 
systems. Our main conclusions are: 

(i) We find no evidence for any anistropy in the satel- 
lite galaxy distribution relative to the major axes of the 
primaries. 

(ii) The projected number density profiles of satel- 



lites brighter than a V-band magnitude of —17 are well 
determined for three separate bins of primary magnitude, 
-21.25,-22,0,-23.0. 

(iii) Apart from the faintest satellites, for which there 
is a slight excess at small galactocentric projected distance, 
all other density profiles are well fitted by projected NFW 
profiles that have been background subtracted to match the 
procedure that has been applied to the data. 

(iv) The concentration of the NFW fits decreases sys- 
tematically with increasing satellite luminosity and is almost 
independent of the luminosity of the primaries (see Fig. [S]). 
Thus, bright satellites have more extended distributions and 
fainter satellites are more centrally concentrated. 

(v) The radial distribution of satellites is dependent on 
the colour and morphology of their primaries. Satellites are 
more numerous around red/early primaries and have more 
extended, lower concentration, distributions (see Fig. (Sjl. 

(vi) The radial distribution of satellites also depends 
on the colour of the satellites. Blue satellites are more nu- 
merous than red satellites at all radii (for the luminosity 
range we probe) and faint red satellites are more centrally 
concentrated (higher NFW concentration) than faint blue 
satellites. Further sub-divided samples show that the con- 
centration of the blue or red satellite profile depends more 
on the colour of the primaries than it does on the colour of 
the satellites. 

As a check of potential systematic effects in our results, 
we have also performed the same analysis using the SDSS 
DR7 dataset. Generally, the results based on DR7 are con- 
sistent with those from DR8, although we do observe some 
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Figure 10. Satellite density profiles split by the colour and luminosity of the satellites. Panels a) and b) compare the profiles of blue 
(blue solid line) and red (red dot-dashed line) satellites for two different bins of satellite luminosity (see legend) . Panels d) and e) show 
the corresponding unsealed profiles as red and blue symbols together with curves depicting NFW fits. The best fitting amplitudes and 
concentrations are given in the legends. For satellites brighter than -20.0, panels c) and f) compare the profiles of red and blue satellites 
around red and blue primaries. The profiles for blue primaries are shown with solid lines and those for red primaries with dashed lines. 
The blue and cyan lines are for blue satellites and red and orange lines for red satellites as detailed in the legends. In panel f) the smooth 
curves show the NFW fits to the corresponding unsealed profiles and again the best fit parameters are listed in the legend. 



differences in the distribution of faint satellites. This is most 
likely due to less accurate photometric reduction and sky- 
background subtraction for DR7 (see Appendix [B] for more 
details) . 

With the advantage of our large and carefully se- 
lected samples, we have discovered a variety of interest- 
ing information about the projected number density pro- 
files, which it has not been possible to quantify clearly 
in previous work. Howeve r, even with a very l imited sam- 
ple, a pioneering study bv lLorrimer et al.l l| 19941 ) found that 
the distribution of satellites is dependent on the morphol- 
ogy of primaries. They found that the number of satel- 
lites around early-type primaries is greater than that about 
late-type primaries and that the concentration of the satel- 
lite distribution is higher around early type primaries. We 
confirm the greater abundance of satellites aro und early- 
type primaries, but contrary to iLorrimer et al. I (fl994i ') we 



find higher concentrations for satellite systems arou n d late - 
type primaries. More r ecently, van den Bosch et al.l (|2005h . 
ISales fc LambasI (|2005l ) and IChen et al.l l|2006l ) studied the 
projected number density profiles of satellites of isolated 
galaxies using larger samples from the Two Degree Field 
Galaxy Redshift Su r vey ( 2dFGRS) and SDSS. Although 
Ivan den Bosch et al.l (|2005l ) cautioned that the profiles from 
2dFGRS were incomplete because of incompleteness in close 
galaxy pairs, their study revea l ed th at satellites followed 
NFW profiles. 'Sale s fc LambasI ^OOE) found that the pro- 
files of satellites depart from a power law at small galac- 
tocentric projected distance, and that they are dependent 
on the colour of the primaries, which is similar to our con- 
clusions (iv) and (v). They also found the distribution of 
satellites to depend on their colour, but argued that this 
may be caused by the correlation between satellites and pri- 
maries. In our study, conclusion (vi) shows that the dis- 
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tribution of satellites not only depends on the properties 
of satellites , but also dep ends on the colou r of primaries. 
IChen et al] |2006, ) and Tollerud et all ()201lh selected sam- 
ples onl y from the SPSS spe ctroscopic catalogue in their 
studies. iToUerud et al.l |20l3) found the 3D number den- 
sity profiles of satellites can be fitted by a power-law with 
a slope p oc r~^'^ . After projecting, the slope of this den- 
sity profile will be close to ours. With a careful treatment 
of interlopers, they fitted the profiles with a power-law form 
and found them to be independent of the luminosity of the 
primarie s. These conclusions are consistent w ith ours. Very 
recently. iLares. Lambas. fc Dommguej (|201ll ) estimated the 
radial density profiles of satellites around primary samples 
brighter than —20.5 and —21.5. They also found the ampli- 
tudes of the profiles depend on the luminosity of the pri- 
maries, and the shapes on their colour. 

The physics of the projected number density profiles of 
satellites involves both the physics of the hierarchical assem- 
bly of dark matter halos and the physics of the galaxy for- 
mation that occurs in these assembling halos. Hence quanti- 
fying these profiles will help constrain both galaxy formation 
models and the nature of the dark matter. We expect that 
our profile results and those of others will be an important 
input into refining theoretical models and the next incarna- 
tion of full N-body/gasdynamic simulation that can resolve 
the physics of the formation of satellite galaxies. 
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APPENDIX A: THE PROJECTED NEW 
PROFILE WITH BACKGROUND 
SUBTRACTION 



The NEW density profile ijNavarro, Frenk. fc Wh"it3 Il996l . 
[l993) is 



p{r) 



5c pc 



(r/r-s)(l + r/rs 



(Al) 



where pc is the critical density, 5c is the characteristic over- 
density of the halo and Vs is a characteristic scale length. 
Conventionally, the scale length is specified in terms of a 
concentration defined as c = r2oo / '"s , where the r2oo radius 
is defined as the radius at which the mean interior density 
is 200pc. With these definitions it follows that 

^3 



5c = 



200 



(A2) 



3 ln(l + c)-c/(l + c)' 

We can integrate along a line of sight to obtain the projected 
surface mass density 

E(i?) = 25c Per.. 



1 



Vr^ - (r/rs) (1 + r/r-,)^ 



dr, 



(A3) 



where R is the projected distance from the centr e of the 
halo. This integral can be solved analytically (jBartelmannl 
Il996h and expressed as 



2ScPcrs 

25cpcrs 
3 

25cpcrs 
(l-x2) 



1 - 



zarctanh 



1 + x 
X < 1, 
a; = 1, 



(A4) 



1 - 



zarctan 



x-1 
1 + x 
X > 1 



where x — R/rs- In our profile measurements, we remove 
the contamination of interlopers by subtracting the mean 
density of galaxies in an outer annulus. This outer annulus 
will also contain genuine satellites that are in the outer an- 
nulus of the density profile. Hence, to compare fairly with 
the measured profiles, we should apply the same background 
subtraction process to the projected NEW profile. We de- 
note the resulting background-subtracted projected NEW 
profile as 

„ 2 /•'2X200 

t{x) = T.{x) - T.(x)xdx. (A5) 

200 J X2QQ 

These background-subtracted profiles are compared to 
their unsubtracted counterparts in Fig. lAll The subtracted 
profiles tend to zero at projected R/r2oo ~ 1.4. 



APPENDIX B: COMPARISON OF RESULTS 
FROM DR8 WITH DR7 

We have performed our measurement of satellite number 
density profiles for both the SDSS DR8 and DR7. This helps 
to quantify the impact on the number density profiles of 
the different sky-subtraction algorithms used to define these 
galaxy catalogues. Looking at images of some of our pri- 
mary galaxies in DR7, there were occasions when close-in 
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Figure Al. A comparison of projected NFW profiles (dasiied 
lines) and background-subtracted projected NFW profiles (solid 
lines) for different values of the concentration, c. The projected 
radius, R, is expressed in units of the r2oo radius, r2oo- 



satellite galaxies existed that were not present in DR8. The 
suggestion is that these are spurious fragments of the pri- 
mary galaxy itself, following inexact subtraction of the back- 
ground sky level. One would expect such a problem to be 
worse for lower luminosity satellites and also if the inner 
radius cut to be considered is reduced beyond our default 
1.5 times the Petrosian Rqo- Fig. IBll shows some illustrative 
results where satellites down to 0.5 times the Petrosian Rqo 
have been included in the profiles around Mc = —23.0 pri- 
maries. There is a tendency for DR7 to have extra low lumi- 
nosity satellites near to the primary, which is not shared by 
DR8. This is particularly evident in the lower panel. Further- 
more, while the DR8 profile is robust to changing the inner 
radius cut, the result for the low luminosity satellite profile 
around bright primaries from DR7 changes significantly. We 
conclude that DR7 contains more spurious fragmentation of 
bright primary galaxies, and that DR8 is preferable for our 
study, both in terms of the reliability of the faint galaxies 
and their improved photometry. 



V 



O 




-1.5 -1.0 -0.5 0.0 
Log(r/0 

Figure Bl. A comparison of profiles based on SDSS DR8 and 
DR7 for a) more luminous satellites, and b) less luminous satel- 
lites. The black (solid) lines are the profiles from DR7. The blue 
(dotted) lines are the profiles from DR8. A cut at 0.5 times the 
Petrosian Rgo radius is used to highlight the difference. 

Besides the physically motivated parameters, we also test 
the parameters of the estimation method. The parameter 
Qp helps us to distinguish genuine satellite galaxies from 
background galaxies by excluding galaxies that are at a sig- 
nificantly different redshift. Panel (d) shows that our results 
are insensitive to reasonable changes in the value of Up. 



APPENDIX C: VALIDATION OF THE 
SATELLITE SEARCH PARAMETERS 

In Fig. ICll we show the effect on the estimated number 
density profiles of varying various satellite search param- 
eters. Panel (a) demonstrates that varying iiinncr between 
0.35 and 0.50 Mpc does not change the profiles signifi- 
cantly, because i?innor ~ 0.35 Mpc is already large enough 
to enclose the whole satellite system for primaries satisfying 
—22.5 < Mp < —21.5. The satellite number density profile 
is similarly robust to changes in 7?outer, which is the outer 
radius for the background region, as shown in panel (b). 
The next panel shows the effect of varying AMfaint , the pa- 
rameter used to determine if a primary is isolated. There 
is a very weak variation of the profile shown in panel (c), 
with primaries allowed to have neighbours with a magni- 
tude difference as small as AMfaint ~ 0.1 having slightly 
more satellites than those with larger magnitude differences 
to their neighbours, more in keeping with the term isolated. 
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Figure CI. The effect on the estimated number density profiles of varying the parameters {/Jinncr, Router, AMfainti "p} from their 
default values, {0.3 Mpc,0.6 Mpc,0.5,2.5}, as indicated in the legends. Some error bars for different datasets have been slightly shifted 
for clarity. 
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